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Abstract 

We present a generally covariant formulation of conformal higher-order viscoelastic 
fluid mechanics with strain allowed to take arbitrarily large values. We give a general 
prescription to determine the dynamics of a relativistic viscoelastic fluid in a way 
consistent with the hypothesis of local thermodynamic equilibrium and the second 
law of thermodynamics. We then elaborately study the transient time scales at which 
the strain almost relaxes and becomes proportional to the gradients of velocity. We 
particularly show that a conformal second-order fluid with all possible parameters in 
the constitutive equations can be obtained without breaking the hypothesis of local 
thermodynamic equilibrium, if the conformal fluid is defined as the long time limit of a 
conformal second-order viscoelastic system. We also discuss how local thermodynamic 
equilibrium could be understood in the context of the fluid/gravity correspondence. 
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1. Introduction 



Viscoelasticity is the property shared by almost all continuum materials, showing elasticity 
at short time scales while behaving as a viscous fluid at long time scales [HE]- In papers 
[H H] , the present authors constructed relativistic viscoelastic fluid mechanics in a generally 
covariant form based on Onsager's linear nonequilibrium thermodynamics. We showed there 
that, for arbitrary parameters in the constitutive equations, viscoelastic fluids thus defined 
behave as standard first-order viscous fluids (obeying the relativistic Navier-Stokes equations 
O [6]) at long time scales [1]. We also showed that the evolution equations are hyperbolic 
for a wide range of parameters due to the elasticity at short time scales. In this sense, a 
relativistic viscoelastic model with such parameters gives a causal completion of standard 
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first-order relativistic fluid mechanics [1]. 

Recently, various models of relativistic fluids with second-order corrections in the deriva- 
tive expansion (called second-order fluid mechanics) have been considered in the analysis of 
heavy-ion collision experiments and also in the study of the holographic duality between the 
long wavelength dynamics of black hole horizons and the dynamics of viscous fluids (see, 
e.g., Ej and references therein). For both cases, fluid systems are well approximated to 
be invariant under conformal transformations. In particular, in [9] and [10] , a model of con- 
formal fluid with full second-order corrections was constructed and shown to be consistent 
with the second law of thermodynamics for a wide range of parameters [TU] . 

A remarkable point for second-order fluids is that their entropy densities generally con- 
tain spatial derivatives of thermodynamic variables. This means that local thermodynamic 
equilibrium is broken for second-order fluids. In fact, in thermodynamics the coordinates 
X = (x°, x) represent a coarse-grained spacetime point, whose temporal and spatial resolu- 
tions we denote by et and eg , respectively [3] . We say that local thermodynamic equilibrium 
is realized at x with resolution (et,es) if a small spatial region around x (of linear size eg) 
at time x° can be well regarded as being in thermodynamic equilibrium at least for a time 
duration of et E This implies that the local entropy in the coarse-grained region around x 
is already maximized for given values of local thermodynamic variables at x (such as the 
energy- momentum density and the charge density), and thus that the entropy density at 
X , s{x), is a function only of these local thermodynamic variables and should not depend 
on their spatial derivatives (which include contributions from nearby material particles). 
Thus, when we consider thermodynamics of second-order fluids, we are generally forced 
to give up the hypothesis of local thermodynamic equilibrium and to use (a variant of) 
extended thermodynamics, where the entropy density can depend on quantities including 
spatial derivatives (such as dissipative currents) [131 EH [IS] . 

In this paper, introducing the elastic strain tensor as one of local thermodynamic vari- 
ables, and utilizing the manifestly Weyl-covariant formulation of conformal second-order 
fluid mechanics developed by Loganayagam [10], we construct conformal higher-order (i.e., 
not first-order) viscoelastic fluid mechanics in the Landau-Lifshitz frame in such a way that 
local thermodynamic equilibrium and the second law of thermodynamics are manifestly re- 
alized. We show that the conformal second-order fluid mechanics of [HI [ID] with all possible 
parameters in the constitutive equations can be fully recovered as the long time limit of 



^We call such spatial regions material particles. 
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our viscoelastic model of second order. Thus, if we define conformal second-order fluid me- 
chanics as the long time limit of a conformal viscoelastic system, conformal second-order 
fluid mechanics can be constructed without violating the hypothesis of local thermodynamic 
equilibrium]^ 

This paper is organized as follows. In section 2, we briefly review a part of the manifestly 
Weyl-covariant formulation of conformal second-order fluid mechanics [10]. In section 3, as- 
suming local thermodynamic equilibrium, we present a general theory describing relativistic 
viscoelastic fluids with large strain. In section 4, we construct conformal higher-order vis- 
coelastic fluid mechanics, based on the manifestly Weyl-covariant formalism. In section 5, we 
investigate the transient time scales at which the strain almost relaxes and becomes propor- 
tional to the gradients of velocity. We there verify our claim that a conformal second-order 
fluid with all possible parameters in the constitutive equations can be obtained as the long 
time limit of a conformal viscoelastic system that satisfies the hypothesis of local thermo- 
dynamic equilibrium. In section 6, we briefly discuss how local thermodynamic equilibrium 
could be understood in the context of the fluid/gravity correspondence, and point out that 
its manifest realization may lead to the viscoelasticity/quantum gravity correspondence. 
Appendix |X] gives the list of the dimensions and the weights of various local thermody- 
namic quantities, and Appendix [B] collects useful formulas in the manifestly Weyl-covariant 
formalism with proofs. 



2. Manifestly Weyl-covariant formalism 

In this section, in order to fix our notation, we give a brief review on the manifestly Weyl- 
covariant formulation of second-order fluid mechanics developed in [TU] (see also [TB| fT7]). 

We consider a ci- dimensional spacetime with background metric g^y of signature (— , +, 

^Wc comment that local thermodynamic equilibrium can also be realized in divergence-type fluid me- 
chanics developed by Geroeh and Lindblom [11] [12] . There, a symmetric traceless tensor ^^i, is introduced 
as a local dynamical variable in addition to the standard dynamical variables ^ = /i/T and = u^^/T (T: 
temperature, u'': velocity, /i: chemical potential). In fact, the causally completed extensions of Eckart's 
fluid mechanics [11] and of Landau-Lifshitz's [12] satisfy local thermodynamic equilibrium in the sense that 
the entropy density depends only on ^ , and ^^^i, . The apparent difference between divergence-type fluid 
mechanics and viscoelastic fluid mechanics is that the additional dynamical variable in the former is related 
to the conserved current (often denoted by A^'^'^), while the strain tensor introduced in the latter has no 
origin as a conserved quantity. A possible relationship between two theories will be commented in section [Sj 
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...,+). A tensor Q^^'" is called a conformal tensor of weight w if it transforms as 

QM.-- = e-0g'M... (2.1) 

under the Weyl transformatiorl^ 

9,u = e-"^ . (2.2) 
We introduce the Weyl connection F^^^ = T^^^ + W^^^, with 

r"^. = \ g'"" [d.gua + d.g,^ - d^g^,) , W\, = g^. A' - 5^ A, - 6^ A, . (2.3) 

-'^(w)^ij^ is invariant under the Weyl transformation ( 12. 2p if the (non-conformal) vector field 
A^ transforms as 

A^ = Al-d^ct>. (2.4) 

Following [To], we introduce Dirac's co-covariant derivative [18] for a conformal tensor Q'^'" 
of weight w as 

= v^Q^^z + 1^1^ g^::: + ■ ■ ■ - q;;: - ■ ■ ■ + ^ a gi:;;; , (2.5) 

which enjoys the following properties [TO] : 

v^Q^,-: = e'''^v'^Q'!^::: if g;;::: = e""^ g (2.6) 

Vxg^^u = , Pa'?''" = . (2.7) 
Note that if a contravariant vector f has weight w = d, the following equality holds: 

V^v^" = V^v^" + {w-d) A^yf" = Vf.v^' (when w = d) . (2.8) 
Similarly, if a (2, 0) tensor Q^'^ is symmetric traceless and has weight w = d + 2 , we have 
V^Q'"' = V^,Q^'' {when w = d + 2 , Q>"' = Q''^' and Q"^ = 0) . (2.9) 

Given the velocity field u'^ for a fluid (having unit weight and being normalized as 
u'^Ufj, = —1), the vector field A^^ can be uniquely determined by requiring the covariant 
derivative of to be transverse {u^VyU^ = 0) and divergenceless {V^u^ = 0) [TO] : 

A" = a" - -^u" , (2.10) 
a — 1 



•^Note that g^jy itself is a conformal tensor of weight w = —2 . The dimensions and the weights of various 
local thermodynamic quantities are listed in Appendix [X] 



4 



where a'^ and are the acceleration and the expansion, respectively: 

a" = VuU^ = u''V^u^' , ^ = VX- (2-11) 

In this paper, we write the (anti-)symmetrization of indices as Q[^y) = (1/2) (Q^,^ + Qvt?} 
and = (1/2) (Q/i!/ — <5i//i)- We further denote the symmetric, transverse, traceless part 
of Q^u by 

Q{,u) = KhlQ^^p) - {h^^Qc^p) V , (2-12) 

where h^^ = g^y + u^u,^ is the metric projected to a surface orthogonal to . One then can 
show [To] that V^u" can be decomposed as 

DX = < + < (2-13) 

withe 

(^tiu = 'Dif.Uu) = V{f,Uy) , = l^lfiUu] = K ^[pUa] ■ (2.14) 

We introduce the Weyl-covariantized Riemann tensor TZ^^x" [10] as the curvature for the 
Weyl connection r^^)^^l2 

Tl^vx" = — S^r^^) ,yx + C?!^r(^) ~ ^w) tip ^(t) i/A + ^(w) lyp '^{w) fiX ^ '^pyX ■ (2.15) 

The explicit form is given by 

Ti^iuXa = RpuXa + J^pu QXa " 4 (J^fif^J^A^f] (^V„^;3 + " ^S'"/?) = S"^*^ "^^^.^Aa • (2-16) 

It is easy to see that the following equality holds for a covariant vector vx of weight w |10] : 

\D^,Vj\vx = Tltiux'' v„ + wJ^^yVx, (2.17) 

J^^. = 5^ A - • (2-18) 
For a contravariant vector of weight Wi, we obtaiii^ 

V,] v'l = -7^^..^ t;;^ + 7^, v\ . (2.19) 



"'Note that '^{y.v) ■ We will use the abbreviation such as {(t^)^i, = cr^ aav and tr((T^) = craO-p^ ■ It 

then holds that (cr")^^ = (cr")(A'i') h^hj^ {a''')(^ai3) for n = 1, 2, . . . , and (w")^,. is h^hf {uj"\af3) for even 
n and h^hf {uj"-)[ap] for odd n. 

^We follow the convention of [TBIIIT] where all the curvature tensors are negative of those in [TU]. Tensors 

of non-calligraphic font are the ones constructed from the affine connection F^^^ (i.e. the ones obtained by 
replacing by V 



P '^j ^ pj 

^Note that wi ^ w + 2 \i ^ g-^^ . 
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We also introduce the following conformal tensors: 



TZ = n^ = R + 2{d- 1) VpA" -{d-l){d- 2) = e^* 7^' 



7111/ 



n 



fii/ 



Gp, + J-^. + (rf - 2) 



^ 3 



pv 



These tensors have the following symmetry properties for their indices [TU] : 



-Ac ;t!/ 

2 t^'^ ' 



Ao- - 



7^ 



(2.20) 
(2.21) 

(2.22) 

(2.23) 
(2.24) 

(2.25) 

(2.26) 
(2.27) 

(2.28) 
(2.29) 

(2.30) 

(1 — 2 " i/*^""^"'^! y""'^ 2{d—l)^"^^J " ^ pvXcF • (2-31) 
Although this is given only in terms of the metric, we can rewrite it as a sum of Weyl- 
covariantized curvature tensors [see Eq. f lB.6p ]: 

4 /_ _R 

l{d 



The Bianchi identity ( [Dx, [D^j, D^] ] + (cyclic) ) Vp = gives [TO 

T^xT^pvp^ + T^pTtuXp^ + T^vTlxp.p^ = , 

Contracting the indices in the first equation, we obtain 



a I'^fj.vX 



fi'Vi/A 



0, 



The last equation is equivalent to 



d-2 



The Weyl tensor is defined as 

4 



C, 



R 



fii/Xa 



C' 



fii/Xa 



n 



4 1^ 



When Ap, is constructed from the conformal velocity field of unit weight as in Eq. f l2.10p 
we have the following formulas [see Eqs. (1B.7P and (]B.2p ]: 

1 



d-2 

u'^TZ^pu^ = -tr(a2) _tr(a;2)_ 



n 



{pv) 



(2.33) 
(2.34) 
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The shear cr^^ and the vorticity u^^ satisfy the following equalities [V^ = w^P^) [see 
Eqs. flRl2D . flRTOj) . dEl, dRT]) and (EID]: 

D^D.a^'^ = 7^<^'^> a^, - ^ J-'^'^a;^, + i p„ 7^ + (rf _ 2) n,) , (2.35) 

P^P^cu^^ = - ^ ^^'^u;^, , (2.36) 

T^u(T^iu = U°'Ca^iul3U^ + ^ _ 2 "^(a"^) ~ ~ ('^^)(mi'> ' (2.37) 

= ^ /i/ -T^Aa - (cTW + wa)^,, . (2.38) 
One can further show that 

= V,V,a^^ - ^ F^^-uj,,, - (rf - 1) V^{F^''u,) , (2.39) 
P^(M,7^'^^/^/) = 2tr(a3) + 6tr(aa;') - 2a'"' u'' C^^.^u^ " ^^''''^/^- 

- ^ J^'^'^.u + V.V.a^'' -{d- 1)P^( J-^^n,) . (2.40) 



3. Relativistic viscoelastic fluids with large strain 

In this paper, we consider a conformal viscoelastic fluid in a rf-dimensional spacetime with 
background metric g^u{x) . We assume in this section that there exists a conserved charge 
(such as particle number) in addition to energy and momentum, and later will ignore the 
charge (or set the corresponding chemical potential to be zero) to simplify discussions and 
expressions. A viscoelastic fluid then has the following local thermodynamic variables [3111]: 

p^(x), n(x), g^,v{.x) , e^,y{x). (3.1) 

Here, p^{x) is the energy-momentum vector, n[x) the charge density, and e{x) = (£:^j,(x)) 
the strain tensor, which is assumed to be spatial {e^^u^ = 0) and symmetric {e^^, = Sy^) 
[3l Hj. We define the proper energy density as 

e(a^) = \/-g^'"(,x)p^{x)py{x) , (3.2) 

which is essentially the sum of the rest mass energy density and the internal energy density. 
Note that are always additive quantities, but e is not in a relativistic theory. 

In this paper, we work in the Landau-Lifshitz frame, where the velocity field u = u'^ is 
defined as = p^/e (note that u'^ = — l)Ill The metric projected to a surface orthogonal 



'^In this paper, we lower (or raise) indices always with g^i, (or with its inverse g^'^)- 
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to is defined as 



hf,u = Ofiu + u^u^ , (3.3) 
and represents the shape of a material [H H] . By writing the velocity field aqf 

u = u>^d, = ^do + ^d, (3.4) 

with the lapse function and the shift functions TV* , the metric g^i, can be expressed with 
the following ADM parametrization: 

ds^ = -N^ [dx'^Y + hij {dx' - N' dx°) {dx^ - dx°) . (3.5) 

Note that h^^ and h'^'^ are written as 

A™ A«A"A" -A..A'\ \ 

\h,o h„j X-h.tN" h„ j \(l (h-^)'j 

where is the inverse matrix of the (d — 1} x [d — 1) matrix h = (hij) . Thus, the local 
volume y/h = ^/de^hij) satisfies the identitjc 

= ^ {h~T ^h., = ^ 5V = ^ h^"" ■ (3-7) 

We assume that local thermodynamic equilibrium is realized at each spacetime point, so 
that the local entropy density s{x) at spacetime point x is given as a function of the values 
of these local thermodynamic variables at the same point x, s{x) = s(p^(x), n{x), g^^{x), 
£fj.u{x))- In order to treat conserved quantities, it is often convenient to multiply densities 
by the local volume y/h , and we write such multiplied densities with placing tilde on the 
original densities. Then, the fundamental relation for the local entropy s = \/h s is given by 

[a ID 



6~s = -^6p,-^6h + ^ rr 6g,, + ^ Si^ 6e; . (3.8) 

Here, T and fi are the temperature and the chemical potential, respectively, and Tq"^ = 
e u^u^+Tq^ is the quasi-conservative energy-momentum tensor with r^^ the quasi-conservative 



^We denote spacetime coordinates by x = (a;^) = (a;°, x') (/i = 0, 1, . . . , 1; i = 1, . . . , d— 1). Note that 
we need not assume u'^ to be hypersurfacc orthogonal when considering only infinitesimal neighbors around 
a spacetime point. 

^To obtain the last equality, we have used the identity h^^" 5{u^Uv) — . 
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stress tensor (tq"^ 



Tq^, Uy = 0)|1j We write the above with (not with e^y) for later 
convenience. S*^ is assumed to be spatial and symmetric (5''"^ = 5**^^) and represents the 
entropic (not energetic) force caused by strain^ If we rewrite and 5g^y in Eq. (13 .Sp in 
terms of 5e , 5u^ and 5h^u by using = eu^ and g^^i, = h^u — Ufi Uy , then terms proportional 
to 5u^ are all canceled, and we obtain another expression for the fundamental relation: 

1 



5~s 



T 



5e 



T 



12 



We assume the densities with tilde to be extensive: 

A5(e, n, h^,, = s{X~e, An, \^/^''~'^ h^,, e';) (VA > 0) . 



(3.9) 



(3.10) 



Then, taking a derivative with respect to A and setting A = 1 afterwards, we obtain the 
following Euler relation, having the same form as that for a fluid: 



e + P — fin 



or 



e + P — fin 
T 



(3.11) 



T 

Here, P = (tr ro)/((i — 1) = g^u Tq'^ / {d — 1) = h^^ r^'^ / {d — 1) corresponds to the pressure of 
an isotropic fluid. From Eqs. ( 13. 7p . (13. 9p and (13. lip , we obtain the fundamental relation for 



1 



T T 2T ° 
together with the Gibbs-Duhem relation. 



TSi^ Ss ,, , 







(3.12) 



(3.13) 



We deflne the order of a local quantity as the order of derivatives necessary to be taken 
in order to construct the quantity from local thermodynamic variables p^, n, g^^, and e^y : 

order 



Pfi ; ; d^u ; ^ fJ-u ; ^ ) y T , ft , Tq , S 



7^ 



liuXrr 1 • • • 





1 

2 



(3.14) 



^^T/^" and Tq'' were denoted by Tj^^^ and r^'^'^ , respectively, in [31ll]. 
We call 5^ the entropic force since it works such as to maximize the local entropy s . This is in contrast 
to the (energetic) elastic force that is exerted on a material particle through the spatial divergence of the 



stress tensor (see Eq. (|4.1ip '). 

^^If we use instead e^i/ as an independent variable, then it should be scaled in the same way as that for 



Ltf 5 ^ Illy 



^2/(^-1), 
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Accordingly, the energy-momentum tensor of a viscoelastic fluid has the following expansion 
in the Landau-Lifshitz frame: 



= eu>'u'' + Tr + rr + r^'' + --- (rf V = ; ^ = 0, 1, 2, ■ ■ ■ ) , (3.15) 
where the subscript 0, 1, 2, ... of a term stands for its order. 

4. Conformal viscoelastic fluids with large strain 
4.1. Definition 

We say that a viscoelastic fluid is conformal if (1) its energy- momentum tensor is traceless 
(T^ = 0) and has weight w = d + 2 and (2) its equations of motion are Weyl covariant. 
We assume that the strain e = (e^,^) has the same weight as that of the metric Qf^^, (i.e., 
w = —2). We further assume that the trace of T'^^ vanishes at each order, obtaining the 
equalities 

e = trro = (rf- 1)P, trr^ = (i = l,2,...). (4.1) 

Thus, the conformal energy-momentum tensor is generically written in the following form: 

T^^ = e u'^u'' + h"'' + T^"''^ 

a — 1 

= e u^^u" + ^ h'^" + T^^"^ + ri""^ + + ■■■ . (4.2) 
To avoid discussions (and expressions) from being unnecessarily messy, we will set /i = in 



what follows 
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For a system of infinite size with free boundary condition, a system which has neither 
energy-momentum transfers nor strains at spatial infinity, a global equilibrium is character- 
ized by the conditions that (i) ds/dp^ = —u^/T be spatially constant and (ii) = [3]. 
The Weyl-covariant expression for the condition (i) is given by 

= /i^" (-] = - P^n^ + j J u'' (in global equilibrium) , (4.3) 



^^For a conformal fluid (in the absence of strain), one can show from the Gibbs-Duhem relation (|3.13p 
that e and n are written in the form e = e T'^/x''"'^ and n — [{d~k)/{d— 1)) eT'^fj,'^^'^^-^ with some constants 
e and k . 
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or 

(i-a) Vf,u,, = a^^ + ujf,^ = . 

(m global equilibrium) . (4.4) 

(i-b) V^T (= h- {d^T + a„ T)) = 

We know from (i-b) that the nonvanishing of spatial covariant derivatives of T indicates 
departure of the system from global equilibrium]^ so that it should be more natural to 
express thermodynamic quantities with the pair (T, h, e) than with (e, h, e) when making 
the derivative expansion around a global equilibrium state. The temperature T is then the 
only dimensionful quantity for a conformal viscoelastic fluid, and we can express various 
thermodynamic quantities as follows [a = (da/j), co = (wq/j)]: 

s = T''-'s{h,e), e = T'e{h,e), S^^ = S^^{h,e) , (4.5) 
= T'^ fi;\h, e) , rr = f^\h, e, a/T, co/T) , ... (4.6) 

From the fundamental relation fl3.12p and the Euler relation (13. lip (with fi = and P = 
e/{d — 1)), we obtain the following equations: 

d _ d ^_ d 



e , 6s = —-6e = - + d del . (4.7) 



rf- 1 ' d-l 2 



4.2. Equations of motion 

The equations of motion comprise (1) the conservation law of T^'^ and (2) the rheology 
equations [20l [2T| [3l H]. The former, V^T^'^ = 0, determines the evolution of the energy- 
momentum p^^ = eu^ (or equivalently, that of T and u^) and describes the motion of material 
particles of a given material, while the latter determine the evolution of the strain e = {s^y) 
and describe the plastic deformation of the material. 

(1) conservation law 

Since T^^ is symmetric traceless and has weight w = d + 2, the conservation law can be 

^■*The acceleration = u^V = u'^ (S^u^ — df^u^) can be written as = N^^ d^N for the ADM 
parametrization p.4p and (|3.5p . Thus, the condition (i-b) can be rewritten as h^da{NT) = 0, which 
is equivalent to di{NT) = 0. Note that To = NT is the temperature conjugate to the energy density 
measured with (not to the proper energy density measured with the local proper time) which includes 
the gravitational potential (Tolman's law), and that it is Tq which becomes spatially constant in global 
equilibrium when the gravitational field exists [31 [TH [H] . 
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(4.8) 



The longitudinal component can be further rewritten as follows: 



(Jj X 

= V^e + T^^""^ a^, . (4.9) 



This determines the time evolution of e . Setting 5 = Vu in the Gibbs-Duhem relation f l3.13p 
[with /i = and P = e/{d — 1)], and using the identity VuQ^y = T^uh^y = , we can also 
write down the equation determining the time evolution of T : 

VuT _lVue d-lS^'^ 
T ~ d e d e 

= - ^ [f'^""^ a^, + id-l) S'^" V^e,,] . (4.10) 

The transverse component of the conservation law can be interpreted in two ways. When 
it is written with the standard covariant derivative , the equation h^JDaT'^^ = gives 
the equation 

= m'^ V,M^ = - - /i^, V,r"^ (4.11) 

e 

which determines the time evolution of and describes how a material particle moves under 
the influence of the stress t'^'^ . The second interpretation is obtained when the equation is 
written in a manifestly Weyl-covariant form: 

= h^"^ V^T""" = h^'^ V^t""" = h^'^ {T'^ f"^) , (4.12) 

or 

= /i^,I)„r°- + rff^^ (/i," ^) , (4.13) 

from which the spatial derivatives of the temperatures, VaT/T , can be expressed in 
terms of other local variables. 



^^See [3] for the derivation of the conservation law V^T^'' = on the basis of Onsager's hnear regression 



theory. 
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Equations fl4.1Up and f l4.13p and the identity V^h^i, = thus allow us to concentrate our 
consideration only on the time evolution of the strain. 

(2) rheology equations 

The rheology equations, on the other hand, have the following form in the derivative 
expansion: 

VuE^, = [V^e^,]o + + [V^e^,]2 + ■■■ . (4.14) 

Their explicit form will be given in section |5] assuming that the strain is small [see Eqs. f lS.lOp 
and f im]) ]. 

4.3. Entropy production and the second law of thermodynamics 

The entropy current consists of the convective part s and the dissipative part s^^-^ , 

si^ = su^ + s[^-^, (4.15) 

from which the total entropy density is defined as 

stot = -u^s^' = s + As {As = -u^ sj'j)) . (4.16) 

Under the hypothesis of local thermodynamic equilibrium (assumed throughout the present 
paper). As is a local function of thermodynamic variables (not depending on their spatial 
derivatives) and can be absorbed into the original entropy density s . We, however, leave 
the possibility of the existence of As in order to make easier the comparison of the results 
to be obtained with those obtained in other references based on extended thermodynamics. 

The entropy production rate is then given by 

V^s^ = V^s^ = Vus + I)^sf^) . (4.17) 

The first term of fl4.17p can be expressed explicitly by using Eqs. fl3.12p and 04. 9 p as 

Vus = ^ P.e + i S^, V^e; = r^^^^ a,, + ^ S^'' V^e,, 

= T"^-! W ^^^^ + s''- v^e^^) , (4.18) 

where we again have used the identity VuQ^u = T^uhfj^u = . The second term of (14.171) can 
be calculated once sj^^-j is given explicitly, whose generic form in the derivative expansion is: 

sjj^ = T^-i (m + n^) (n^n^ = 0) 

= T^-^ [ (mo + mi + ma + ■ ■ ■ ) + ^0 + ^1 + ^2 + ■ ■ ■ ] • (4-19) 
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Here, Uq = since a dimensionless transverse vector cannot be constructed only from h^^, , 
e^iy and . We will set rfiQ = in the following discussions since rfiQ can be absorbed into 
s = s/T'^~^ . We leave the possibility of the existence of ffii {i = 1,2, . . .) for the reason 
stated below Eq. IKT6\i . 



5. Conformal viscoelastic fluids at transient time scales 
5.1. Small strain expansion 

We denote the shear (i.e. traceless) component of the strain tensor e = (efj,^) by cs = {^f^u) = 
{£(f_iu)) , and decompose the strain as 



^M^ = ^M^ + ^^(tr^) V- (5-1) 



The rheology equations then take the following form in the derivative expansioni^ 

V^e^^ = const. T + const, + ■ • ■ , (5.2) 
Vutie = const. T tr 5 + const. T tr(£:|) + const. tT{e^a) + const. T"^ tr((T^) H . (5.3) 



As time elapses, the strain gets relaxed and the rate of change becomes small. If the time 
scale of observation is sufficiently long, then the left hand sides of Eqs. (15. 2p and (15. 3 p become 
negligible compared to the right hand sides, and we obtain^ 



const. T + ■ ■ • , (5.4) 
tr e — J- const. tr(£:|) + const. tr^es cr) + const. tr((T^) + ■ • • 

~ const. tr(cr2) H . (5.5) 

Thus, in order to investigate a conformal viscoelastic fluid at the transient time scales where 
T^uS^y is small but still not negligible, it is convenient to expand the equations of motion 



^^An explicit paranietrization is given in Eqs. (|5.10p and (|5.1ip . 

^^Note that, unlike the shear part ()5.4p . the trace part, tre, cannot contain terms proportional to the 
gradient of velocity, = V^^u^ , in the derivative expansion, since "d is not a conformal scalar. This is 
consistent with the absence of bulk viscosity in conformal fluids. 
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based on the following semi-long time (SLT) order: 





SLT-order 


e, T, . . . 

C/^i/) ^fiv^ ^^u^ ■ ■ ■ 

tie, TZ(f^u), ■ ■ ■ 



1 
2 



(5.6) 



For the rest of this section, we investigate a conformal viscoelastic fluid at the transient time 
scales to second SLT-order in the constitutive equations. 

We first expand the energy-momentum tensor 
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d 



(5.7) 



to second SLT-order as 

We also expand the entropic force caused by the strain, S^'^ = T'^ S^'^ , as 

1 



(5.^ 



and the rheology equations a^l^ 



tr£ + ■■■)/^'^^ (5.9) 



ci a^^ + C2T-i(a2)<'"^> +c3T-i(cra;)<^^> +C4T-i(u;2)(/.^> 

+ rfi T e^^ + d2T {el) + 4 (^5 a) ^'^'^^ + d^ [es uj)^^"^ + ■■■ , 
VutTe = foT tie + fiT tr(£|) + /a tr (£5 ^) + /s tr (a^) + /4 T^^ tr (a;^ 
+ hTtiiel) + /e tiielcr) + /^T"! tr(55t^2) 



(5.10) 



+ fsT-' tT{escru) + fgT-^ tries 



fioT 



-2 tiia'] 



-2 tr^aa;2^ 



+ /12 M° C^^^,)p u^ + /i3 a^'^ m° 

+ /i4 7^(^,) + /i5 7^(^,) 

+ /16 + /i7 + ■ ■ • . 



(5.11) 



^^As can be seen from Eqs. p.33p and p.37p . only two are independent among the symmetric, transverse, 
traceless tensors of weight two; Vucr^,, , u°''Ra(tiu)p'^^ ^ '^{p.'^) ' '^"^aij.yfiU^ We will take 'R(^v) and 

as such two in the following discussions. 
^^The equations generalize those given in [20l[2Tl|4j (where only ci and di_ are nonvanishing) . 
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Then, by using Eq. f l4.18p . VuS can be expressed explicitly as foUowsl^ 

VuS = T'^[ii di tr(4) + {-bi + ii ci) T-^ tries cr) - ai T'^ tT{a^) 
+ {£2 di + h d2) tr(4) + (-62 + £2 ci + £1 4) T-^ ti{el a) 
+ (-63 + £1 C2) tr(e5 (T^) + (&4 + ^1 C3) tr(£5 a w) + £1 C4 T^^ tr(£5 c^^) 

- as tr(a3) - 04 T^'^ ti{a u^) + £1 C5 

- a5 T-3 a'^'^ M° + £1 C6 T'^ 7^(^.) - ag T'^ a^'^ 7^(^,) ] . (5.12) 

As for the dissipative part s^^^ in the entropy current = su^ + s^^^ , we need to expand 
it to second SLT-order|f^ 

.f,) = T'^-i { [AT"i ti{esa) + (A2/2) T^^ ti{a') + (A3/2)T"2 tr(u;2) 
+ (A4/2) tr(e2 a) + tr(55 a^) + T"^ tr(55 a 

+ tr(£5C^') + ^8T-2 7^ + ■■■]M^ 

+ V^e'^s'' + Aio P.a'^'^ + A^ T"' P.w^" + A12 T'^ M,7^"''/^/ + ■■■}. 

(5.13) 

As mentioned before, the scalar in front of represents the correction As to the original 
entropy density s ; As = Stot ^ s. Local thermodynamic equilibrium is inevitably broken 
when any of these coefficients Ai {i = 1, . . . , 8) do not vanish. Although we eventually set 
= (i = 1, . . . , 8) later, we leave them for a while for comparison with other references. 

The entropy production rate can then be written in the following form [note that P^s'^ = 
V^sf" due to Eq. fl^ ]: 

P^s^ = (s + sji)) = VuS + D^sj'd) 

= T'^ { £1 di tr(e|) + (-61 + £1 ci + Ai di) T"^ tr(e5 a) + (-ai + A^ ci) T'^ ti{a^) 

+ (£2 di + li d2) tr(4) + (-^2 + ^2 ci + Ai d2 + ^4^1 + li 4) T"^ tr(e| cr) 
+ (-63 + A4 ci + £1 C2 + c/i + Ai d-i - Ai) T-^ ii{es a^) 

+ (&4 + ^1 C3 + A di - ^4) ir{es + (£1 C4 + - Ai) T'^ tr(e5 a;^) 
+ (-as + C2 + ci - + 2A12) tr((T^) 

^°Since (es)"^ = eg , cr'^ = ct , w'^ = — cj , we have tr(e5(Tw) = (1/2) tr((e5(T + CTes)^) = (1/2) tr(e5((7aj — 
wcr)) . 

^^We have neglected terms proportional to T'^-^ tre or T''^^ tr(e|) since they can be absorbed into 
s [see a comment below Eq. (|4.19p ] . We also have neglected terms proportional to F^^" Uy since it is at least 
of third-order derivative [see Eq. (jB.lSP ], 
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+ (-04 + AtCi + Ai Ci- A2- 2A3 + 6A12) T-^ tr(a 
Ai 



+ (£1 C5 + Ai) + (-05 + A,c, + A2- 2A^,) T~' a^'' C^^^^^p 



+ Kl C6 + 



c/-2 



° d-2 
{2As + A,o + Au)T-^V,V,a^^'' 



(5.14) 



Here we have used Eq. (12.401) and neglected terms proportional to 'D^^{J^^yU'^) or 'D^jV^uj 
= —{{d — 3)/2^J^fj_^u'^^ since they are at least forth-order derivatives [see Eqs. (IB.ISP and 
( IB.17P ]. We also have used the fact that T can be treated as being covariantly constant to 
this order (P^T ~ 0) since both VuT and h^V^T are at least of second SLT-order [see 
Eqs. fl410|) and ( 1413|) ]. 

Now, if we set 



Ao = 



2A8 + Aio + A 



112 



0, 



the entropy production rate takes the following bilinear form in our approximation: 



D^s^ = T'^ V^l M V^"" 



(5.15) 



(5.16) 



with 



M 



2 



2 

-ai + Ai ci 



*8x8 / 



at"' /T 

4) 



(4)(A'-> 



(cr2)(M>')/T2 
(aa;)<'"'>/T2 



(5.17) 



The second law of thermodynamics is then certified if the coefficient matrix M. is positive 
semi-definite. 

The conformal fiuid mechanics of [10] is obtained by setting the parameters as follows: 



ai = r]i, 02 = 

05 = ^72 + ^76 , «6 
d-A 



V2 + V4., 03 = -2 ?73 , a4 = -?72 + ?75 

V2 



A2 
A12 



V2 



d-2 



Vg 



d-2 

A3 = 

Ai = Ai 



bi Cj di fi 







V5 



d + 2 



V6 



2 V"^ 'd-2 

^5 = ^6 = ^7 = Ag 



A: 



1,2,---), 



' 2(rf-2) 
Aio = 0. 



A 



11 



^5 + 3?76 
2(d-3) ' 

(5.1^ 
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A calculation based on the fluid/gravity correspondence shows that A2 , A-^ , Ag 
P, Uni US] , and thus local thermodynamic equilibrium is broken for such conformal fluids. 
In the next subsection, we show that the breakdown of local thermodynamic equilibrium 
can be avoided if a conformal fluid is always defined as the long time limit of a conformal 
viscoelastic system. 



5.2. Long time limit and second-order fluid mechanics 

We now consider a conformal higher-order viscoelastic system when the time scale of obser- 
vation is much longer than the relaxation times of the strain: 



s 1 



(5.19) 



As in first-order viscoelastic systems (see, e.g., our second-order viscoelastic system 
comes to behave as a second-order viscous fluid (which is conformal now). In fact, the 
rheology equations f IS.lOp and fIS.lip gives the following equations in the long time limit 
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^ - — \ciT-'a>''' + C2T 



-2 (0-2^^1.) 



C3T ^{cru 



C5 



+ d2 {el) ^^"^ + dsT~' {es (r) + d.T^^ {eg co^^"^ + ■ ■ ■ ] 



di 

di 
1 



C2 ^ _ dsci 

di di \diJ d\ 



Vdi di J 



-2 (auj)^^"'^ 



2 I 2 



di di 



tie ^ --[h tr(4) + ti{es<y) + /sT"^ ^(a^) + f^T^^ ti{uj'' 

Jo 



■] 



r/i 









fo di fo 



rp~2 ^ /„2^ 



o- 



^T"2 tT{uj^) + . 
Jo 



(5.20) 



(5.21) 



Substituting this into the constitutive equations (15. 7p and (15. 8p . we obtain the energy- 
momentum tensor in the long time limit: 

^(Cng) = «o T'^ {9^" + d u'^u'') + 7 T'^"' 



, 02 C?? + bi (Ci rfi - ^2 - C2 rf? ) + 62 di - 63 Ci ._2 , 2n (ui.) 

+ T (a)^'^) 

Q3 c^i + ^1 (ci 4 - C3 C?i) - 64 Ci cil 



04 - C4 ^^_2 ('^2>i(mi'> 



^^Note that the long time hmit of shear strain, Eq. (|5.20p . has the same form as the additional dynamical 
variable in divergence- type conformal fluid mechanics (see Eq. (98) of [Hj). 



18 



+ "^^^7^^"^ C-^^^^f' u, + 7 T'^"^ 7^<^'^> + ■ ■ • . (5.22) 

d\ di 

Here, we have set e = (rf — 1) (^o- constant) since e becomes constant in the long time 
hmit (Stefan-Boltzmann law). Equation fl5.22p has the same form as the energy- momentum 
tensor of a generic conformal second-order fluid jTOl UHl EI] • 

Interestingly, even if we start from a viscoelastic system with manifest local thermody- 
namic equilibrium [i.e. a system with Ai = {i = 1, ... ,8) so that As = 0], all terms can 
appear in Tj^^^^^ with any possible values, in a way consistent with the second law of ther- 
modynamics. In fact, there are no constraints on parameters in the constitutive equations 
in order for the strain e^i, to be converted to spatial derivatives (such as a^i, and u^y) in 
the long time limit [see Eqs. f l5.20p and fl5.2ip ]. Furthermore, we can also understand the 
appearance of spatial derivatives in the entropy density of a conformal fluid as a result of 
the same conversion mechanism, which is now applied to the entropy density with manifest 
local thermodynamic equilibrium: 

s = s(p^, Qf,^, e^y) = s{p^, g^y, 0) + const. T'^"^ tr(£:|) + const. T'^^^ tre H , (5.23) 

that transmutes in the long time limit into the one with spatial derivatives: 

S(iong) = s(p^, 0) + const. T"^-^ tr(a2) + const. T'^-^ tr(a;^) + const. T'^'^ 7^ + ■ ■ ■ . 

(5.24) 

Thus, even though the hypothesis of local thermodynamic equilibrium holds at short time 
scales for our viscoelastic system, it is seemingly broken when the system is observed at long 
time scales and is treated as a viscous fluid. 



6. Conclusion and discussions 

In this paper, we defined conformal higher-order viscoelastic fluid mechanics. We wrote down 
the equations of motion in such a way that the evolution is consistent with the second law of 
thermodynamics. We further showed that any conformal second-order fluid with arbitrary 
parameters in the constitutive equations can be obtained by taking the long time limit 
of a viscoelastic conformal fluid, without violating the hypothesis of local thermodynamic 
equilibrium. 

On the other hand, if one trusts the fluid/gravity correspondence, the entropy current 
of a conformal fluid can be computed in the gravity side. The result [lOl [16] shows that the 
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total entropy density s^ot contains spatial derivative terms with nonvanishing coefficients, 
and thus we know that local thermodynamic equilibrium is violated even at short distance 
scales for such conformal fluids that have gravity dualso 

As was argued in [3], even when local thermodynamic equilibrium is realized for a system 
with resolution [e^, eg), spatial derivative terms are naturally induced in the total entropy 
density (as the entropy functional in the language of [3] or as in Eq. (15.241) ) if we observe the 
system at larger scales in both the temporal and the spatial directions. Since viscoelastic 
fluids allow a description with manifest local thermodynamic equilibrium, we expect that 
the fluid/gravity correspondence is an already coarse-grained correspondence between vis- 
coelastic fluid mechanics and a more microscopic description of gravity. If this is the case, it 
then should give an important clue to finding fundamental degrees of freedom in quantum 
gravity to try to formulate such "viscoelasticity/quantum gravity correspondence." A study 
along this line is now in progress and will be reported elsewhere [22] ■ 

As another direction of future research, it would be interesting to apply viscoelastic fluid 
mechanics to the phenomenology of heavy- ion collision experiments. In fact, relativistic 
viscoelastic model gives a causal completion of relativistic fluid mechanics (the latter being 
defined as the long time limit of the former) [1], and thus it is tempting to assume that 
there is a phase of viscoelasticity prior to the stage of viscous fluidity. Then, it is important 
to investigate how elasticity at short time scales affects the dynamics of states right after 
collisions. In particular, one should investigate whether elasticity drives the system to an 
ideal fluid more rapidly than in the standard second-order fluid mechanics, as has been 
observed in divergence-type fluid mechanics [12]. 
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A. Weights of local thermodynamic variables 

The Weyl weight of a {p, q) tensor Qtl-.-v^ of dimension A is given bytf = A+ p — g. We 
list below the dimensions and the weights of various local thermodynamic quantities. 





U-liiiCliblUli -LA 


Wclg,liL IxJ 


(Ji LLci 
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AT / A 

N/A 







-1 
i 
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— z 
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tie 
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-2 













-2 
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-2 


2 
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-2 


2 






2 





2 




T — d A — d A 


2 





2 






2 


2 


2 






d 


d- 1 










d 


d 








S s(^6, £ ^u) 


d- 1 


d- 1 








T = (ds/de)-^ 


1 


1 








s^' = su^" + sf^^ 


d- 1 


d 


N/A 


N/A 




d 


d + 2 


N/A 


N/A 



(A.l) 



B. Useful formulas 



In this appendix, we prove a few useful formulas which are used in the main text. 
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For the Weyl-covariantized Riemann tensor fl2.15l) - fl2.18l) . the following equality holds: 



- {<y''')tiu - {<y (^)t,u - <y)txu - {^'^)tJiv - '^^^{a^y + u^^) , (B. i) 



where h^^, = g^^ + u^^Uy and = . By decomposing Eq. (IB.ip into the trace part, the 



symmetric traceless part, and the antisymmetric part, we obtain the following identities!^ 

11^13 = - tr((T2) _ tr(a;2) , (B.2) 

T^a{i,u)p = (0-^)(m!-> + i^'^){f^u) + T^u(^i,u , (B.3) 

i hP hj Tp„ = [auj ^uja)^y^ V^oJ^v , (B.4) 

where we have used the relation IZ^aiup) = ^^laQais and u" 7la[p.u]i3 = (1/2) -^^i/ [see 
Eqs. fl2:23|) and fl2:25|) ]. 

The Weyl tensor (I23TD . 

4 _ R_ 

2{d-l] 

can be rewritten as a sum of Weyl-covariantized curvature tensors with the use of Eqs. fl2.16p . 
(Km and (Km -. 

This exhibits that C^^ao- is a conformal tensor of weight —2. Since the tensor L^^^^ = 
4<^f„^/.][.5^] satisfies ^"L^^^^^^m/^ = -5f^5^^ , we havfl 

Capyp {=u°' C„ (^^) ^ u'^ ) =u°"Ra (^^) /3 u'^ - 7^(^^) . (B.7) 



^^We can show that Eqs. (|B.2p - (jB.4p are equivalent to the well-known evolution equations of a^y, and 
LOfj^i, (such as the Raychaudhuri equation) by using the following equations: 

^2 



d-1 



^^Note that Ca^^/^ u*^ = ^^j,) ^ since Capi./3 = C^i^^a and u^u^^Ca^^p = = Ca^.^/^u'^u^ . 
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Multiplying Eqs. (IB.3P and (IB.4p by a^'^ and u^'^, respectively, and using Eq. (IB.7p . we 
obtain the following formulas: 

tr(fr^) + tr(a u') + a^" , (B.8) 

1 
2 

We can also show 



, , .2\ 



T^^'uj^,, = -2 tr(a a;^) + uj>^- v^co^, . (B.9) 



V.V^co'^'' = ^ [V„ V,] u^'' = - 7^[^.] u;^^ + ^ J"^, u;'^'^ 

= -^^^'^^M^, (B.IO) 

= n^^^^a,, - '^F^'^u,,, (B.ll) 

d 3 

1 Q 



1 „ „ . . , „n . 



= n^^^^a^, + i I)„7^ - ^ ^^^a;^, + (d - 2) [j"^, u^] . (B.12) 

We can show that the spatial vector J-'^jj, u'^ vanishes up to third-order derivatives!^ 

J'^.u u'' = - [a^, + u^J) a" - h^, Vua" + ^ h; = 0{e') . (B.13) 

In fact, from Eqs. (14 .Qp and (14. lip , we find 

-Vue = Wu\^e + -^'d = 0{e^), a'^ = - 3 /i'^" 9„ In e + ©(e^) , (B.14) 
e a — 1 a 

and thus have 



K. Vua" = - ^ V V. (/i"" 5, In e) + 0{e' 



- ] h^u iVuh"'') dalne-] /i " Vuda In e + 0(6^) 

- I V (Vn(w"w")) In e - ^ V V,V, In e + ^ /i/ (V,m^) 9^ In e + 0{e^) 
^ V 9,^9 - + O a, + 0{e') . (B.15) 



^^We denote terms of n"^ and higher order derivatives by 0(e") . 
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Similarly, using Eq. (IB.Mp . we can show that the spatial components of J-'^y vanish up 
to third-order derivatives: 

1 "d 

= -k h;hf V[„ [h^] d. In e) - u:,. + 0(6^) 



from which we have 



- 2 ^i^^- V„ In e - ^/^- + 0{e') 

0{e') , (B.16) 



^^^a;^, = 0(e^). (B.17) 
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